Generic Finite Size Enhancement of Pairing 
in Mesoscopic Fermi Systems 
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The finite size dependent enhancement of pairing in mesoscopic Fermi systems is studied under the assump- 
tion that the BCS approach is valid and that the two body force is size independent. Different systems are 
investigated such as superconducting metallic grains and films as well as atomic nuclei. It is shown that the 
finite size enhancement of pairing in these systems is in part due to the presence of a surface which accounts 
quite well for the data of nuclei and explains a good fraction of the enhancement in Al grains. 



Since long it is a well-known fact that in certain finite size 
Fermi systems the gap is increased substantially from its bulk 
value. Such systems are, for instance, ultra small supercon- 
ducting metallic grains, of great present actuality [ij, yl |3[, 
and thin films |4j, |5|, la, |7fl but also superfluid atomic nuclei 
1 8, 9]. There have been theoretical studies in the past on 
the size dependence of pairing in the abovementioned sys- 
tems I Hi ll ll 1 1 2L 1 1 311 . To our knowledge for the condensed 
matter systems no satisfying explanation has been found fT3l 
whereas for the nuclear systems large scale Hartree-Fock- 
Bogolioubov (HFB) calculations for nuclei have recently 
somewhat clarified the situation flUl . 

In this investigation we will set a rather limiting frame: we 
assume that BCS theory is valid and that the pairing force 
v(r) is size independent. These are, of course, very severe 
restrictions, and obviously, other size dependent features may 
be present in reality. We will consider simplified systems: 
First we study metallic grains and films in a hard wall poten- 
tial using the standard schematic constant matrix element ap- 
proximation with an adjustable strength parameter and a cut- 
off given by the Debye frequency. It will be shown that this 
model acounts for a good fraction of the experimental size de- 
pendence. Second we apply the previously developed pocket 
formula to the mass number dependence of nuclear gaps. We 
will see that our simple theory describes the mass number (^4) 
dependence of nuclear pairing quite well. In all cases only the 
spin singlet channel shall be considered. 

Let us first present our general approach. As already men- 
tioned, we want to base our consideration on the validity of 
BCS theory. In finite systems the gap equation can there- 
fore be written in the standard form where the states 
\n) are the eigenvectors of the single particle Hamiltonian 
h = p 2 /2m* + V (r) with V (r) the (phenomenological) 
single particle potential and m* = m* (r) the effective mass: 



A„ = - 



in n v n n 



A n >/2E n 



(1) 



eigenvalues of h, i.e., h\n) — e n \n), the pairing matrix 
element (nn\v\n'n') contains the time reversed states |n), 
and the chemical potential /i for finite systems is determined 
by the "particle number (AO condition": N = ^2 n |(l — 
(e„ — fj,)/E n ). This model, though quite schematic, will al- 
low us to develop the essential features of the size depen- 
dence of pairing. One further important hypothesis, as al- 
ready mentioned, is that the pairing force from which the ma- 
trix elements in Q are constructed, does itself not depend 
on the size of the system. Still the matrix elements, via the 
wavefunctions, will be size-dependent. One guesses that the 
other important sources of mass number dependence in Q 
are the single particle spectrum, respectively the level density 
g (e) = J2 n ~ £ »)' an d tne cnem i ca l potential ijl. 

We, at first, will apply a statistical approach ltl4[l5ll . This 
essentially consists of replacing the single particle density ma- 
trix \n)(n\ by its value averaged over the energy shell 1 15]. 



Pe„ = -n)]% - e n >)\n')(n' 



S(e n - h). 



(2) 

An asymptotic expression for p £ri can then be derived 
using the semi-classical method by Balian-Bloch for infinite 
hard wall potentials 111611 or the Thomas-Fermi (TF) or equiva- 
lently Strutinsky averaging method for smooth potentials 
Recognising that the two body wavefunctions (fi^rm) in 
the pairing matrix elements can be written as (rir2\rin) = 
(ri\n)(n\r2), we can pass to the continuum limit and write 
forQ 



A(e) = - / de'g(e')v(e, e')A{e')/2E{e'). (3) 
The averaged pairing matrix element is given by 

v(e, e')=[[ dYoT'ff'v (p - p') S(R - R'). (4) 



In Q E n are the quasi-particle energies with E\ = 
(e„ — fjb) + and the single particles energies e n are the 



where dF — dRdp / (2ir) 3 and v(p) is the Fourier transform 
of the pairing force, / = f e ( R, p) is the Wigner transform 1 9 ] 
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of p £ in (|2ji and a prime on V and / means that all variables 
should be replaced by primed ones. The size dependence of 
the gap parameter A = A (e = fi) is then contained in the 
corrections to the bulk values of g(e), v(e, e') , and /i. 

Let us first evaluate A for the case of metallic grains and 
films. The electrons be confined by an infinite hard wall poten- 
tial of arbitrary shape. As usual in condensed matter physics, 
we approximate the attractive electron-electron interaction by 
a delta function pseudo-potential with a cut-off in energy sym- 
metrically on both sides of the Fermi energy \l of the or- 
der of the Debye frequency u>d- In the bulk the pairing ma- 
trix element is therefore given by (k — fc|v|fc' — k') = ^J 2 - 
for |efe — fi\, \ekf — fi\ < u>p> and zero otherwise and V is 
the volume of the system. For a finite size grain our main 
task will be to evaluate the pairing matrix elements @ for 
this case. The expression of the level density g(e) in terms 
of volume, surface, and curvature contributions is well known 
since long fl^l . For the matrix elements we also will em- 
ploy the Balian-Bloch method 1 16] using the method of im- 
ages. To lowest order the distribution functions in @ are 
given by f e (R, p) oc 8 (e — h 2 p 2 /2m) which is the bulk ex- 
pression. In order to obtain the correction term, we transform 
back into coordinate representation, f e (R, p) — > p t (r, r'), 
and then replace z' by —z', the z-direction being the one per- 
pendicular to the surface. Back into phase space one obtains, 
f e (R,p) = g(e)- 1 [6 (e - h 2 p 2 /2m) + Sf] with 



Sf = -5( Pz ) 



2m/h 2 



cos(2Rk e (p x ,p y )) (5) 



where k e (p x ,p y ) = (ff (e - ^ (p 2 x + 2 .Since 

/<= (R,p) is normalized to unity, one obtains from (|5}, in in- 
tegrating over phase space, the classical result for the level 

density g (e) = ^ (f?) 1 yfW - ^ (^)El- An impor- 
tant point to be realised is that the volume V and surface S 
correspond to the borders of the hard wall. Since the density 
is diffuse at the surface, the relevant matter volume Vm < V 
is therefore given by the wall delimitation which encloses the 
correct number of particles. The relations between V, S and 
V m ! S m are worked out in fl7il and are to lowest order given 



by V=V 



M ' 



hd + ■■■ and S=Sm 
at the Fermi energy then becomes : 



The level density 



9f = g(e = fi) 



V M 2m 
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+ ...). (6) 



We remark that the sign of the surface term is now positive, 
that is, for a given volume V m the level density is enhanced 
by the presence of a diffuse surface which, in fact, is the usual 
situation. With (|5jl and the definition of 17(e) it is, in consid- 
ering that (Sf) 2 also contributes to order straightforward 
to evaluate the pairing matrix-element @. In the case of our 
delta force, its Fourier transform is a constant and one obtains, 



v(e,e)= — 



^ 7r min(k ei k e i) S 
+ 4 V 
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n min(k e , k e ') S m 37r Sm 
4 k t k t , V M Sk F V M 



(7) 



We therefore see that, contrary to the level density, the ma- 
trix element vp — v(/i, fi) diminishes in absolute size in 
the presence of a surface. All ingredients are now prepared 
and one can solve the gap equation Q for instance numeri- 
cally. However, there exists a well known and accurate an- 
alytical solution which is more interesting in. The result 
is A = 2lud exp ( v ^ gp ^ . Inserting g F from (|6jl and Vp 
from @ into the above expression, we notice that the prod- 
uct Vpgp does not depend on the surface. However, one also 
has to account for the compression effect due to the surface 
tension which increases the chemical potential or respectively 
the Fermi momentum, and thus gp. Finally this leads to an 
enhancement of the gap for low system sizes. Elaborating one 
obtains k F = k F (l + f pry-JH. where fcj? stands for the 
bulk value. Inserting into the expression for the gap one ob- 



tains 



A = A 



(8) 



where Vp and As stand for bulk values. One clearly sees 
that the gap becomes enhancend as the size of the system de- 
creases. 

It is fortunate that formula (|8j can be tested on a very early 
quantum mechanical solution of JO for a slab lllill . In this 



case one has 



V M 



— 4- where L is the film thickness. In 



1H2I1 the constants in were chosen —v^gp = 0.3 and 
kp = 0.84 x WPcmT 1 . It can be seen from Fig.l that our 
pocket formula passes on average well through the quantum 
mechanical values fT3l . 




FIG. 1: Dependence of the gap, for the case of a superconduct- 
ing homogeneous film, on the film thickness L. The saw tooth line 
corresponds to a quantum mechanical calculation 1 12], whereas the 
smooth curve corresponds to formula (8). The horizontal line rep- 
resents the bulk value As for Aluminium. The dots represent the 
center of gravity of the triangles in which they are lying (a crude way 
to estimate an average of the quantal results) 
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n w 2R[A] 2R[A] A[K] A[K] eq.jj [K] 



60 


41.49 


40.83 


0.00 


0.00 


1.34 


80 


45.73 


45.06 


0.00 


0.00 


1.31 


100 


49.30 


48.64 


0.00 


0.83 


1.28 


200 


62.22 


61.55 


0.95 


1.18 


1.22 


300 


71.26 


70.60 


1.00 


1.18 


1.19 


400 


78.46 


77.79 


1.00 


1.16 


1.17 


500 


84.53 


83.86 


1.00 


1.15 


1.15 


1000 


106.54 


105.87 


1.00 


1.12 


1.12 



TABLE I: Number of levels in the window (n w ), size (2R), (2R) 
and gap (A), (A) without and with surface correction, respectively. 
The gap obtained using eq.^8} is also given. 



gap is larger than the one we have calculated and correspond- 
ingly R c is smaller. However, a precise estimate of the effect 
is difficult. The gap can also be calculated from the exact so- 
lution of the picket fence model (see |21]). It turns out that 
this 'quantal' definition of the gaps yields, around the phase 
transition region, substantially larger values than those from 
the mean field BCS theory, again enhancing the ratio A/ As- 
The quantal values of A also can be obtained for sizes quite a 
bit smaller than R = R c of BCS theory. We therefore think to 
have isolated an important enhancement mechanism of pair- 
ing in metallic nano-grains, stemming from the presence of 
a surface. Other effects, like e.g. the size dependence of 
the phonon spectrum, should be taken into account to obtain 
quantitative agreement with experimental data. 



In refs. [3, IH 0, it is indicated that in the case of 
Al grains one obtains, with respect to the bulk, an enhance- 
ment for the critical temperature T c by roughly a factor of 
two for a grain diameter of 45A. For a spherical grain with 
Vm = 4 ^ one obtains J^- = However, grains are 
rather pancake shaped than spherical |2, 1911 . For an oblate 
ellipsoid with short diameter half the one of a sphere with 
the same volume the increase of Sm/Vm is 44 %. Proba- 
bly grains are even triaxial (see ref.|19|, fig. 2) and we take 
Sm/Vm = 9/(2R) which corresponds to a 50 % increase 
over the spherical case. Taking in (|8j the bulk values for Al 
that is kp = 1.75A -1 and -wf =0.168, we obtain from 
© for A/A B an enhancement — 30 % at 2R — 45A which 
is a sizeable fraction of the experimental value. However, in 
such small grains the electron levels are discrete and it is well 
known |9] that the gap equation has no solution, if the av- 
erage level distance d ^> A^- We therefore solved the gap 
equation Q for the picked fence model (equally spaced levels 
with Kramers degeneracy) JJJ] for cud = 395 K which is the 
value for Al. The number of levels nw in the window 2u>d 
was estimated to be i) n w = 2ojDgp if we take only the low- 
est order term in l|6) and ii) nw — 2uDgF when including 
the surface correction to the level density (and the one coming 
from fi, see above). For the dimensionless interaction constant 

we take —A = v F g F = vfgf (l + fprf^-), with .gf as 
above. In this way we also can calculate A/Ag quantally 
in the picket fence model. We find that A/Ag raises from 
A/A B = 1 for R = oo to A/A B - 1.2 at 2R ~ 60A, fol- 
lowing quite accurately our pocket formula. For smaller grain 
sizes the solution of the gap equation quickly breaks down, the 
critical size occurring at 2R C ~ 40A. The situation is sum- 
marized in Table 1. It therefore seems within our schematic 
model that one can only reach a moderate enhancement of 
20% - 30% depending on whether or not one believes into a 
continuation of the increase into the pair-fluctuating regime. 
Several comments are, however, in order: Equal level spacing 
is the most unfavorable situation which can exist . Usually 
a certain percentage of grains have some symmetries which 
can enhance the gap (see ref. 1 20]). Therefore on average the 




N ( number of neutrons ) 



FIG. 2: Average nuclear gaps as a function of neutron number N 
along the valley of /3-stability of the nuclear chart. The experimental 
points have been taken from 1 10]. Broken line: the asymptotic value 
A_b=0.37 MeV to which the full line converges. 

In nuclear physics it is well known since decades that pair- 
ing is stronger in lighter nuclei than in heavier ones. An em- 
pirical formula A = 12/y/A with A = N + Z the sum 
of neutron (N) and proton (Z) numbers had been used in 
the past to fit the data |8, |9J|. However more recently Sat- 
ula et al. 1 22] pointed out that the data used so far to ex- 
tract the gap values were overestimated and contaminated by 
the Jahn-Teller effect l22ll . A new analysis using the filter 
A = | \E£ +1 + Eq -1 - 2Eg] for neutron number N odd 
only, Eq being the measured binding energies of nuclei, re- 
vealed that the mass number dependence of A is substantially 
weaker than the Yl/\J~A— law. In nuclear physics it is com- 
mon use to solve the gap equation (1) either, as for the metal- 
lic grains, also using a <5-force pseudo-potential with a cut- 
off |H 01 or more sophisticated finite range forces are em- 
ployed for the matrix elements in (1) not necessitating any 
cut-off. One of the best tested and successful forces of the lat- 
ter type is the Gogny D 1 S force 12311 . In principle for nuclei 
it is more appropriate to work with smooth potentials like the 
Woods Saxon or harmonic oscillator potentials and to use for 
the average density matrix on the energy shell (3) the well 
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known Wigner-Kirkwood ^-expansion |9]. This procedure 
is, however, more cumbersome and does not lead to such a 
handy formula as (|SJl. For space reason we cannot present 
this here and it will be published separately in the future. 
For the time being we will also use (|8j for finite nuclei as a 
generic formula. In nuclear physics the convention is such 
that -Vp = ^9- = G and gp 1 = where the level den- 

sity parameter a = ^- ^? b ~a A MeV -1 . An average value 

n k F 

from Skyrme and Gogny forces is a ~ A. MeV -1 . Atypi- 
cal value for G which can be found in the literature | 9| 12411 i s 
G = f MeV. We also checked, using the methods of 1 15], 
that this latter value is compatible with the Gogny D1S force 

On average nuclei are spherical and then = where 

R = tqA^ is the nuclear radius. The product kfro = (^g-)i 
is a universal number and then, besides Ab, all constants in 
l|8} are fixed also for the nuclear case. The bulk value of 
the gap is a quantity which in nuclear physics is quite un- 
certain because the mass number range of nuclei is too small 
to extrapolate to infinite nuclear matter without the guidance 
of a reliable formula. We expect l|8} to be such an expres- 
sion which allows to pin down Ag within certain limits. In 
Fig. 2 n we show that a good fit to the data with the above 
values for a and G is obtained with Ab = -37 MeV. Us- 
ing for a — y^MeV -1 which is obtained with m = m* and 
which is the standard Fermi gas value used in phenomeno- 
logical models, the fit yields Ab = 0.45. This gives a 
slightly flatter but still acceptable curve than the one shown 
in Fig. 2 and shows that formula (|8jl, for the nuclear case, is 
quite robust. These values for A b are of the same order of 
magnitude as the asymptotic value Ab — 0.58 MeV calcu- 
lated from the D1S force Izill . In Fig. 2 the ^4-dependence 
has been converted into an Y-dependence via the relation 
A — N = A/ (1.98 + 0.0155.4!) which defines the valley 
of stability of the nuclear chart 1 25] . Therefore for nuclei the 
pocket formula gives a very satisfying reproduction of the 
data and we thus conclude that it contains the essentials of the 
physics. 

In conclusion, we isolated in this work an important and 
generic enhancement factor of pairing in finite Fermi sys- 
tems. This stems from the surface corrections to their respec- 
tive bulk values of level density, pairing matrix element, and 
chemical potential. We derived a pocket formula for the en- 
hancement factor A/ Ab which is very general and depends 
exponentially on the ratio surface to voulume of systems of 
arbitrary shape. It remains valid for level spacings d < 1.4A 
because for larger spacings the solution of the gap equation 
breaks down. Our theory explains satisfactorily the average 
experimental mass number dependence of nuclei. For Al 
grains we obtain within the picket fence model a maximun en- 
hancement of A /A b ~ 1 -2 at a grain diameter of ~ 6 nm. We 
checked that the situation is similar for the case of Sn grains 
fl^l . This estimate is based on BCS theory. We, however, ar- 
gue that in a more realistic theory the corresponding gap may 
exist for smaller grains yielding a more important fraction of 



the experimental results. Other effects mentioned above can 
give additional enhancements. Studies in this direction are 
planned for the future. 
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